Block magnetism coupled with local distortion in the iron-based spin-ladder BaFe2Se3 
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Magnetism in the insulating BaFe2Se3 was examined through susceptibility, specific heat, resistiv- 
ity and neutron diffraction measurements. After formation of a short-range magnetic correlation, a 
long-range ordering was observed below Tn 255 K. The transition is obscured by bulk properties. 
Magnetic moments (|| a) are arranged to form a Fe4 ferromagnetic unit, and each Fe4 stacks anti- 
ferromagnetically. This block magnetism is of the third type among magnetic structures of ferrous 
materials. The magnetic ordering drives unusually large distortion via magnetoelastic coupling. 
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PACS numbers: 72.20.-i, 75.25.-j, 75.30.Cr, 75.40.-s, 75.80. -|-q 



Since its discovery, research on iron-based supercon- 
ductivity (SC) has become the main stream in condensed 
matter physics^. The interplay between structure, mag- 
netism and SC is one of most intriguing subjects of 
this field. The common theme is on the 2D plane of 
iron, which has been realized in the ZrCuSiAs (1111), 
ThCrsSia (122), anti-PbO (11) and CuzSb (111) struc- 
tures. To gain further insight into the mechanism of SC, 
realization of, for example, different crystal/magnetic 
structures is desired. A recent prominent example is the 
245 system^, in which a novel type of magnetic structure 
in an insulator has been discovered. Regarding crystal 
structure, achievement of SC over various dimensions is 
also ambitious. As spin ladders in copper oxides shed a 
new light on the mechanism of SC, a study on an ana- 
logue with ladder geometry among ferrous compounds is 
highly interesting. Research on that compound can pro- 
vide unique opportunity to study electronic correlations 
and importance of dimensionality for SC. 

BaFe2Se3 is the first realization of ladder geometry 
among iron-based compounds^. Unlike most of parent 
compounds of the iron-based superconductors, BaFe2Se3 
is an insulator down to the lowest measured temperature. 
The structure is novel, comprising of FcSe4 tetrahedra 
with channels which host Ba atoms [Fig. 1(a)]. Four- 
fold coordinated Fc^"*" ions extend along the 6-axis and 
forms two-leg spin-ladder structure within the 6c-plane. 

Here we report study of the first iron-based spin-ladder 
compound BaFe2Se3 based on high-quality polycrys- 
talline sample. High-resolution powder neutron diffrac- 
tion reveals a 3D magnetic order below Tn ^ 255 K. 
However, bulk properties such as magnetic susceptibility, 
electrical resistivity and specific heat show no anomaly 
at the transition. Magnetic moments (~ 2.8 /ib at 5 K) 
are arranged to form a FC4 ferromagnetic (FM) unit per- 



pendicular to the ladder direction, and each Fe4 stacks 
antifcrromagnetically. This block magnetism, as for in 
recently discovered 245 systems^, is the realization of the 
third type of magnetic structures among iron-based ma- 
terials. Rietveld refinements on neutron diffraction data 
also reveal unusually large local lattice distortion driven 
by the magnetic order. 

High-quality polycrystalline BaFe2Se3 was grown by 
the solid state reaction. We first synthesized the precur- 
sor BaSe3. Stoichiometric amounts of elemental Ba shots 
and Se powder in an alumina crucible were sealed into an 
evacuated quartz ampoule. It was slowly heated up to 
700 °C, kept for 12 hrs, and then quenched to room tem- 
perature. The obtained BaSe3 and Fe with a 1:2 molar 
ratio were mixed together and put into an alumina cru- 
cible inside a quartz tube. It was reacted at 700 °C for 96 
hrs with an intermediate regrinding. The powder x-ray 
diffraction using Cu-Ka radiation indicates no trace of 
impurity phases, however, it turned out that there is a 




FIG. 1: (color online), (a) Crystal and (b) magnetic struc- 
tures of BaFe2Se3 at T = 5 K. Cuboids with dashed lines 
indicate crystallographic unit cells. 



TABLE I: Atomic positions within Pnma of BaFe2Se3 atT — 
300 K determined by Rietveld analysis {x^ = 2.3). Lattice 
constants are a = 11.93178(15) A, 6 = 5.43747(6) A, and 
c = 9.16476(12) A. Isotropic Debye-Waller factor ((7iso) is 
employed. 



Atom 


Site 


X 


y 


z 


t/iso (A^) 


Ba 


4c 


0.1862(5) 


1/4 


0.5154(10) 


0.051(5) 


Fe 


8d 


0.4941(3) 


0.0000(6) 


0.3525(3) 


0.044(2) 


Sel 


4c 


0.3595(4) 


1/4 


0.2300(5) 


0.059(5) 


Se2 


4c 


0.6283(3) 


1/4 


0.4923(6) 


0.034(4) 


Se3 


4c 


0.3997(4) 


1/4 


0.8155(5) 


0.061(5) 



tiny amount of FerSeg impurity by magnetization mea- 
surements. The electrical resistivity was measured by the 
standard four-probe method. The magnetic data were 
collected by using a SQUID magnetometer. The spe- 
cific heat was measured by utilizing a commercial set-up 
(Quantum Design, PPMS), where the relaxation method 
is used. Neutron powder diffraction data were collected 
on high-resolution ECHIDNA diffractometer at ANSTO 
and WAND at ORNL, with A = 2.4395 and 1.4827 A, re- 
spectively. Diffraction patterns were obtained between 
T — 5 and 350 K in a closed-cycle-refrigerator. We 
employed group theoretical analysis to identify magnetic 
structures that are allowed by symmetry. 

Table |T] summarizes the structural parameters at 
T = 300 K determined by Rietveld refinement on pow- 
der neutron diffraction data. Consistent with other 
reports^"—, the crystal structure is accounted for by or- 
thorhombic space group Pnma. By keeping site occu- 
pancy variable for each atoms in the refinement, the 
result indicates stoichiometric ratio within error bars, 
Bao.96(4)Fe2.oi(2)Se3 02(7)- Temperature dependences of 
lattice parameters are provided in Fig. 2(a), showing 
systematic change throughout temperature regime. 

Electrical resistivity p shows an insulating behavior 
[Fig. 2(b)]. Note that p exhibits no anomaly in the 
whole temperature range measured. A fitting with an 
activation-type formula p = Aexp(A/T) yields a small 
activation energy A = 0.13 eV, comparable with Rei^. 
More precisely, p deviates from the activation-type be- 
havior, instead it shows a variable-range-hopping (VRH) 



type behavior. The data is fit to p oc exp(i3 -I- (^)T+d)j 
where d denotes the dimension of the VRH model. The 
VRH in ID model yields better fit rather than higher 
dimensions [inset to Fig. 2(b)]. This indicates that the 
compound is classified as a Mott- Anderson insulator, and 
is in the vicinity of metal-insulator (MI) transition. To 
further examine the conduction behavior, measurements 
on a single crystal are now under way. 

Accurate magnetization measurement under small ap- 
plied field is limited by a tiny amount of FM impurity, 
FctScs {Tq ^ 460 K^). We therefore deduce the suscepti- 
bility X = M/H under p,qH = 5 T and subtract the one 
at 1 T to minimize effects from FeyScg. The key issue 
on the susceptibility is that it shows a weak temperature 
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FIG. 2: (color online). Temperature dependences of (a) lat- 
tice parameters, (b) electrical resistivity (p), (c) magnetic sus- 
ceptibility (x) and its inverse, and (d) specific heat (Cp). The 
insets of (b,c,d) show a fit with ID variable-range-hopping 
model (see text), enlargement of bifurcation, anomaly around 
Tn (left top) and Cp jT against with a polynomial fit (right 
bottom), respectively. 



variation in the paramagnetic phase (above Tn) and does 
not obey the simple Curie- Weiss law [Fig. 2(c)]. One pos- 
sible explanation for this is that a simple localized spin 
picture in a Mott- Anderson insulating state is invalid in 
the present system, which is close to the MI transition as 
anticipated by the small charge gap. Another possibility 
could be that the unquenched orbital angular momentum 
is present in Fe^+ ions owing to the prominence of the 
spin-orbit coupling in an orbital degenerated state. In 
contrast to the earlier report^ indicating an emergence 
of weak ferromagnetism below Tn, our high- field suscep- 
tibility data show no prominent anomaly over the whole 
temperature range. This possibility is also excluded by 
our powder neutron data. As inset to Fig. 2(c) enlarged, 
X in zero-field-cooled (ZFC) and field-cooled (EC) se- 
quences show small hysteresis. This could refiect glassy 
behavior caused by small randomness, but the origin of 
it is still unclear. From the data of p [Fig. 2(b)] and x 
[Fig. 2(c)], we exclude the possibility of SC behavior of 
this compound, distinguished from Refi^. 

Specific heat Cp shows no anomaly in the whole tem- 
perature measured [Fig. 2(d)]. Given no peak formation 
at around Tn [left-top inset to Fig. 2(d)], entropy re- 
lease associated with the magnetic ordering turned out 
to be surprisingly small. A possible reason for this is 
that entropy release gradually occurs in a wide temper- 
ature range due to large fiuctuations characterized in a 
low-dimensional system as will be discussed later. Cp/T 
versus T^ plot at low temperatures indicates a negligi- 
ble small T-linear contribution, as well as a small de- 
viation from Cp oc behavior [right-bottom inset to 
Fig. 2(d)]. The lack of 7T term is in agreement with 
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FIG. 3: (color online), (a) High-resolution neutron powder 
diffraction pattern of BaFe2Se3 at 5 K obtained on ECHIDNA 
with Rietveld refinement (solid lines). The calculated posi- 
tions of nuclear and magnetic reflections are indicated (green 
ticks). The bottom line gives the difference between observed 
and calculated intensities. The inset shows temperature evo- 
lution of a magnetic reflection with short- (with black curve 
fit) and long-range correlations. 



the Mott insulating nature, and inadequate fit with only 
I3T^ implies an existence of soft phonons. The data can 
be well reproduced by a function Cp ~ f3T^ + with 
/? = 1.89 mJ/mol and 5 = 0.017 mJ/mol K^. If 
it is assumed that all the components of the /3 term 
come from the phonon contribution, using the formula 
/3 = 12tt^ NR/bel, {N being the number of atoms in a for- 
mula unit, 6), the Debye temperature (0d) is estimated 
to be 183 K. This value is considerably small and further 
supports the existence of soft phonons. 

To gain microscopic information of magnetism, we per- 
formed high-resolution powder neutron diffraction. The 
obtained data is well fit based only on BaFe2Sc3, and 
FerSeg impurity has so small volume fraction (< 0.3 %) 
that the secondary phase is not taken into account in the 
refinements. Additional resolution-limited peaks associ- 
ated with magnetic ordering were observed below 250 K 
data. Figure 3 shows a diffraction pattern taken at 5 K 
along with the Rietveld refinement. To identify the mag- 
netic structure, we applied representation analysis. All 
the magnetic peak positions can be accounted for by a 
propagation wave vector, (1/2,1/2,1/2). Basis vectors 
(BVs) of the irreducible representations (irreps) for the 
wave vector (Table ini were obtained using the SARA/i 
code^. There arc 24 BVs in total, belonging to distinct 
two irreps Fi and F2. Each irrep has 12 BVs, in which 
only 4 of 8 Fe atoms in a unit cell are allowed to possess 
magnetic moments, and every BV describes the relation 
of moment direction of 2 atoms either parallel or antipar- 
allel along one crystallographic axis. Both irreps in this 
case are required to participate as a corepresentation to 
let all 8 atoms have magnetic moments. In a certain mo- 
ment direction along either the a, b or c-axis, there are 
2 choices for 4 combinations, namely 16 patterns. We 
assumed that the moment is parallel or antiparallel to 



TABLE II: Basis vectors (BVs) of irreducible representations 
(irreps) for the space group Pnma with the magnetic wave 
vector Qm = (1/2, 1/2, 1/2). Superscripts show the moment 
direction. Columns for positions represent #1: {x,y,z), #2: 
{x + 1/2, -y + 3/2, -z + 1/2), #3: {-x +l,y- 1/2, -z + 1), 
#4: {-x + l/2,-y + l,z + l/2),#5: {-x + 1, -y + 1, -z + 1), 
#6: (-X + 1/2, y- 1/2,^4-1/2), #7: (:e, -j/ + 3/2, 2), and 
#8: {x+l/2,y,-z + l/2). 
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one axis and has the same coefficient of 4 BVs, and sort 
out all 48 patterns to determine the magnetic structure 
by comparing i?-factors. The best fit with i?,„ag = 8.54 
% is a combination of BVs, -04, -^lo, V'le and ?/'22j which 
is consistent with Ref^. The estimated moment at 5 K is 
2.75(2) /iB per an Fe^"*" site. The 5 K magnetic structure 
corresponding to the best fit is shown in Fig. 1(b). Mag- 
netic moments (|| a) are arranged to form a Fe4 FM unit 
in the 6c-plane, and each Fe4 stacks antiferromagneti- 
cally. For inter-ladder coupling on the ac-plane, block 
spins that are AF correlated in one neighboring ladder 
and ferromagnetically correlated in the perpendicular di- 
rection. The next nearest-neighbor (nn) ladders interact 
antifcrromagnetically. The determined magnetic struc- 
ture, described as "block magnetism" , is reminiscent of 
one observed in insulating ^2Fc4Sc5 (A — K, Rb, Cs, 
Tl)2, with vacancy-ordered 2D Fe plane. 

Despite diverse structural and physical properties, 
there have been so far reported only several types of mag- 
netic structures among iron-based materials. The com- 
mon theme of magnetism and SC is on a square lattice of 
Fe with either pnictogen or chalcogen. In cuprates SC, 
parent compounds as Mott insulators show a simple AF 
pattern, so-called checkerboard magnetic structure. This 
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FIG. 4: (color online), (a) A ladder geometry in BaFe2Se3. 
Solid lines surround unit cells. Temperature dependence of 
(b) estimated magnetic moment of the LRO, (c) ratio between 
lengths along the ladder {u, v) and leg (w) directions, and (d) 
representative angles (ai,Qf2) of Se-Fe-Se. 



is favored by strong AF nn interactions, yet has not been 
observed in any ferrous compounds. The most common 
spin structure among the ferrous is single-stripe magnetic 
structure. This features AF nn spins in one direction and 
ferromagnetically correlated in the perpendicular direc- 
tion. Stripes extend along nn bonds as favored by next 
nn AF exchange interactions. The single-stripe is ob- 
served in the parent compounds of the 1111, 122 and 111 
materials. Double-stripe is also widely observed in the 
11 systemaiS and insulating La202Fc20Se2^, in which 
stripes extending again along nn bonds in a structure 
that is favored by AF third nn interactions. The mag- 
net CaFe4As3 with similar FeAs strips but connected in 
3D has an incommensurate variant of the nearly double- 
stripe structurei^ii^. Finally, block magnetism has been 
recognized since the discovery of the 245 systems, in 
which nn 4 Fe atoms connected FM and Fe plane is 
occupied by that block interacting AF with a separa- 
tion of Fe-vacancy. Note that this structure has mag- 
netic moments perpendicular to the Fe plane, whereas 
the other structures have within plane moments. There 
is a close relation between the magnetic structure and 
moment size; single-stripe typically has typically ^ 1 ^b, 
double-stripe > 2 and block magnetism around 3 /ib- 
This classification is also the case for BaFc2Se3, which 
has large 2.8 ^b) magnetic moments perpendicular to 
the ladder-plane, belonging to the block magnetism. 

Temperature dependence of the determined magnetic 
moment is shown in Fig. 4(b). Together with ECHIDNA 
and WAND data, Tn is roughly estimated to be 254.6 ± 
2.4 K by an order parameter-type fit. However, previous 
Mossbauer experiment^ reports no anomaly at around 
Tn, instead clarifies the appearance of magnetic hypcr- 



finc splitting below 100 K. This distinct behavior could 
originate from the different timescale of the experimen- 
tal technique; neutron typically has 10~^^ sec timescale, 
being faster than that of Mossbauer (10"'' sec). Such a 
gradual slowing down of spin fluctuations is sometimes 
observed in low-dimensional magnets. If entropy release 
between 100 K and Tn was sufficient, the obscurity of 
it at the LRO might be uncontradicted. Moreover, we 
observed a broad magnetic peak above Tn at slightly off 
position (Q ~ 0.68 A'^) from 0.73 A'^ of the LRO [in- 
set to Fig. 3]. This diffusive peak with the estimated 
correlation length of 19(5) A, could also reflect the low- 
dimensionality of the material. Elucidation of spin dy- 
namics and short-range correlation realized by BaFe2Se3 
would be important, and will be done shortly. 

Despite the systematic lattice shrinkage as tempera- 
ture decreased [Fig. 2(a)], significant local lattice dis- 
tortion driven by the magnetic order is observed. Fig- 
ure 4(a) depicts one ladder, where the summation of 
lengths, w -|- w is equivalent to the lattice constant b, and 
w is the width of the ladder. Above Tn, u/w and v/w 
is completely overlapping each other and takes slightly 
off value from the exact ladder (1.000) as shown in Fig. 
4(c). They start to deviate after the magnetic order oc- 
curs, v/w shows a small increasing or constant behavior, 
whereas u/w is surely decreasing on cooling. This is at- 
tributed to change of the y position of Fe. We note that 
the atomic displacement along the 6-axis of Fe below Tn 
is as large as approximately 0.01 A, which exceeds 10~^ 
A reported in ordinary materialsi^. This is in contrast 
to the 245 systems which show a shrinkage to the cen- 
ter of Fe4 direction as favored by nn FM interactions^^. 
The distortion in BaFe2Se3 is also evidenced by change 
of Se-Fe-Se bond angles [Fig. 4(d)]. Even above Tn, 
an FeSe4 tetrahedron is distorted, and the magnetic or- 
dering enhances it by elongating the lattice along the 
c. This splits e orbitals and stabilizes — y"^ against 
3^2 _ j,2 ^ ^jjgj-g tjjg direction of z orbital is perpendicular 
to the ladder plane. BaFe2As2 shows also similar temper- 
ature evolution of As-Fc-As angles across tetragonal to 
orthorhombic structural transitiortii. However, the an- 
gle change in BaFe2As2 is half as large as in BaFe2Se3. 
Given these facts of significant lattice distortion inher- 
ent, the — y^ orbital becomes firmly stabilized be- 
low Tn. Again, no symmetry change occurs at all in 
the temperature range from 5 to 350 K in BaFe2Se3. 
The orthorhombic BaFe2Se3 undergoes such a strong 
change due to magnetoelastic coupling without break- 
ing the high-temperature symmetry. Instead it gives rise 
to the unusually large atomic displacement of Fe within 
the unit cell, and it can be consistent with soft-phonon 
mode inferred from the specific heat results. It has been 
reported the orthorhombicity ratio in BaFe2As2 is up to 
3 X 10^"^, and it is systematically suppressed by electron 
doping to finally reach to SC— . BaFe2Se3, on the other 
hand, has almost twice large distortion. Since this type 
of distortion is not preferred by SC in iron-based materi- 
als, applying pressure and doping studies are now under 
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way to accomplish SC. 

In summary we have examined Fe magnetism in the 
spin- ladder geometry realized in BaFe2Se3. A magnetic 
short-range correlation is detected, followed by a 3D long- 
range ordering. The transition is, however, obscured by 
bulk properties measurements. The determined magnetic 
structure is so-called block magnetism, where || a mo- 
ments are arranged to form a Fe4 FM unit, and each 



Fe4 stacks antifcrromagnetically. The magnetic ordering 
also drives unusually large lattice distortion, which SC in 
ferrous compounds does not prefer. 

We thank M. Isobe for technical assistance and K. 
Okazaki for discussion. The work is supported in part 
by Grant-in- Aids for Scientific Research (2340097) and 
the Japanese Society for Neutron Science. 



^ For a review, see e.g., G.R. Stewart, Rev. Mod. Phys. 83, 
1589 (2011). 

^ F. Ye, S. Chi, W. Bao, X.F. Wang, J.J. Ying, X.H. Chen, 
H.D. Wang, C.H. Dong, and M. Fang , Phys. Rev. Lett. 
107, 137003 (2011). 

^ H.Y. Hong and H. Steinfink, J. Solid State Chem. 5, 93 
(1972). 

* A. Krzton-Maziopa, E. Pomjakushina, V. Pomjakushin, D. 
Sheptyakov, D. Chernysliov, V. Svitlyk, and K. Conder, J. 
Phys.: Cond. Mat. 23, 402201 (2011). 

^ J.M. Caron, J.R. Neilson, D.C. Miller, A. Llobet, and T.M. 

McQueen, Phys. Rev. B 84, 180409(R) (2011). 
^ B. Saparov, S. Calder, B. Sipos, H. Cao, S. Chi, D.J. Singh, 

A.D. Christiansen, M.D. Lumsden, and A.S. Sefat, Phys. 

Rev. B 84, 245132 (2011). 

H. Lei, H. Ryu, A.L Frenkel, and C. Petrovic, Phys. Rev. 
B 84, 214511 (2011). 

* A. Okazaki, J. Phys. Soc. Jpn. 16, 1162 (1961). 

^ A.S. Wills, PhysicaB 276-278, 680 (2000); program avail- 
able from www.ccpl4.ac.uk. 

° W. Bao, Y. Qiu, Q. Huang, M.A. Green, P. Zajdel, M.R. 
Fitzsimmons, M. Zhernenkov, S. Chang, M. Fang, B. Qian, 
E.K. Vehstedt, J. Yang, H.M. Pham, L. Spinu, and Z.Q. 



Mao, Phys. Rev. Lett. 102, 247001 (2009). 

D.C. Free and J.S.O. Evans, Phys. Rev. B 81, 214433 

(2010) . 

P. Manuel, L.C. Chapon, LS. Todorov, D.Y. Chung, J.- 
P. Castellan, S. Rosenkranz, R. Osborn, P. Toledano, and 
M.G. Kanatzidis, Phys. Rev. B 81, 184402 (2010). 
Y. Nambu, L.L. Zhao, E. Morosan, K. Kim, G. Kotliar, P. 
Zajdel, M.A. Green, W. Ratcliff, J. A. Rodriguez-Rivera, 
and C. Broholm, Phys. Rev. Lett. 106, 037201 (2011). 
M. Reissner, W. Steiner, and H. Boiler, Hyper. Int. C 5, 
197 (2002). 

S. Lee, A. Pirogov, M. Kang, K.-H. Jang, M. Yonemura, T. 
Kamiyama, S.-W. Cheong, F. Gozzo, N. Shin, H. Kimura, 
Y. Noda, and J.-G. Park, Nature 451, 805 (2008). 
W. Bao, private communication. 

M. Rotter, M. Tegel, D. Johrendt, I. Schellenberg, W. Her- 
mes, and R. Pottgen, Phys. Rev. B 78, 020503(R) (2008). 
M.G. Kim, R.M. Fernandes, A. Kreyssig, J.W. Kim, A. 
Thaler, S.L. Budfko, P.C. Canfield, R.J. McQueeney J. 
Schmahan, and A.I. Goldman, Phys. Rev. B 83, 134522 

(2011) . 



